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Abstract 

Two families of (i-orthogonal polynomials related to su(2) are identified and studied. The 
algebraic setting allows for their full characterization (explicit expressions, recurrence relations, 
difference equations, generating functions, etc.). In the limit where su(2) contracts to the 
Heisenberg-Weyl algebra f)i, the polynomials tend to the standard Meixner polynomials and 
d-Charlier polynomials, respectively. 

Keywords: Generalized hypergeometric polynomials, d-orthogonal polynomials, su(2) algebra, spe- 
cial functions. 

1 Introduction 

We identify in this paper two families of (i-orthogonal polynomials that are associated to 51i(2). 
When available, algebraic models for orthogonal polynomials provide a cogent framework for the 
characterization of these special functions. They also point to the likelihood of seeing the corre- 
sponding polynomials occur in the description of physical systems whose symmetry generators form 
the algebra in question. 

d-orthogonal polynomials generalize the standard orthogonal polynomials in that they obey 
higher recurrence relations. They will be defined below and have been seen to possess various appli- 
cations |6[[7l [13]. Recently, two of us have uncovered the connection between (i-Charlier polynomials 
and the Heisenberg algebra. Here we pursue this exploration of d-orthogonal polynomials related 
to Lie algebras by considering the case of su(2). Remarkably, two hypergeometric families of such 
polynomials will be identified and characterized. 

1.1 d-Orthogonal polynomials 

The monic d-orthogonal polynomials P n (k) of degree n (P n (k) = k n + . . .) can be defined by the 
recurrence relation |10| 

d 

P n+1 (k) = kP n {k) - a n ,n-„Pn-^k), (1.1) 
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of order d + 1 with complex coefficients a njm ; the initial conditions are P n = if n < and Pq = 1. 
It is assumed that a n ^ n -d 7^ (non-degeneracy condition). 

When d = 1, it is known that under the condition c n 7^ 0, the polynomials satisfying three-term 
recurrence relations 

P n+ i(k) = kP n {k) - b n P n {k) - c n P n ^(k), 

are orthogonal with respect to a linear functional a 

(a,P n (k)P m (k)) = 5 nm , 

defined on the space of all polynomials. 

When d > 1, the polynomials P n (k) obey vector orthogonality relations. This means that there 
exists a set of d linear functionals cij for i = 0, . . . , d — 1 such that the following relations hold: 

(o-j, P„P m ) =0, if m > dn + i, 

(<7i, P n Pdn+i) t^O, if n ^ 0. 



1.2 d-Orthogonal polynomials as generalized hypergeometric functions 

Of particular interest are d-orthogonal polynomials that can be expressed in terms of generalized 
hypergeometric functions (see for instance [HO [5]). These functions are denoted p F q and are defined 

by 



7? 



{a p } 1 
{b q Y c 



where (m)o = 1 and (rn)i~ = (m)(m + 1) • • • (m + k — 1) stands for the Pochhammer symbol. In 
the case where one of the a^s is a negative integer, say a\ = —n for n S N, the series truncates at 
/i = n and we can write 



1+sFq 



-n, {a s } 1 
{b q } ''c 



/i=0 



(1.2) 



If one of the bi's is also a negative integer, the corresponding sequence of polynomials is finite. 

The classification of (i-orthogonal polynomials that have a hypergeometric representation of the 
form (|1.2p has been studied recently in [3]. The results are as follows. 

Let s ^ 1 and let {a s } = {a±, . . . , a s } be a set of s polynomials of degree one in the variable k. 
This set is called s-separable if there is a polynomial ir(y) such that Y\i = i(cLi(k)+y) = [QiLi a i(k)} + 
7r(y); an example of such s-separable set is {k e « , x = . . . , s — 1}. If the set {a s } is s-separable, 
then there exists only 2(d + 1) classes of (i-orthogonal polynomials of type (|1.2p corresponding to 
the cases: 

1. s = 0, . . . , d — 1 and q = d; 

2. s = d and g = 0, . . . , d — 1; 

3. s = q = d and 1; 

4. s = 9 = d + 1 and c = 1 and £S Oi(0) - ^=1 ^ ^ N - 

Examples of d-orthogonal polynomials belonging to this classification have been found in U Q3] . 
We shall here provide more examples that are of particular interest as well as cases that fall outside 
the scope of the classification given above. 
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1.3 Purpose and outline 

We investigate in this paper d-orthogonal polynomials associated to su(2). We shall consider two 
operators S and Q, each defined as the product of exponentials in the Lie algebra elements. We 
shall hence determine the action of these operators on the canonical (N+ l)-dimensional irreducible 
representation spaces of su(2). In both cases, the corresponding matrix elements will be found to be 
expressible in terms of <i-orthogonal polynomials, some of them belonging to the above-mentioned 
classification. The connection with the Lie algebra su(2) will be used to fully characterize the two 
families of (i-orthogonal polynomials. The limit as N — > oo, where su(2) contracts to f)i, shall 
also be studied. In this limit, the polynomials are shown to tend on the one hand to the standard 
Meixner polynomials and on the other hand to (i-Charlier polynomials. 

The outline of the paper is as follows. In section 2, we recall, for reference, basic facts about the 
su(2) algebra and its representations. We also define a set of su(2)-coherent states and summarize 
how su(2) contracts to the Heisenberg-Weyl algebra in the limit as N — > oo. We then define the 
operators S and Q that shall be studied along with their matrix elements. The biorthogonality 
and recurrence relations of these matrix elements are made explicit from algebraic considerations. 
Results obtained in |14j and [15] concerning the Meixner and d-Charlier polynomials shall also be 
recalled. In section 3, the polynomials arising from the operator S are completely characterized. The 
result involve two families of polynomials for which the generating functions, difference equations, 
ladder operators, etc. are explicitly provided. The contraction limit is examined in all these instances 
and shown to correspond systematically to the characterization of the Meixner polynomials. In 
section 4, the same program is carried out for the operator Q; the contraction in this case leads to 
(i-Charlier polynomials. We conclude with an outlook. 

2 The algebra su(2), matrix elements and orthogonal polynomials 

In this section, we establish notations and definitions that shall be needed throughout the paper. 
2.1 su(2) essentials 

2.1.1 The su(2) algebra and its irreducible representations 

The Lie algebra su(2) is generated by three operators Jq, J + and J_ that obey the commutation 
relations 

[J+,J_] = 2J , [J ,J±]=±J ± . (2.1) 

The irreducible unitary representations of su(2) are of degree N + 1, with N 6 N. In these repre- 
sentations, Jq = Jo and J± = J=p, where f refers to the hermitian conjugate. We shall denote the 
orthonormal basis vectors by \N,n), for n = 0, . . . , N. The action of the generators on those basis 
vectors is 

J+\N,n) = V \n + Tj(N - n) \N,n + 1) , J_ \N, n) = ^n(N -n+l)\N,n - 1), (2.2) 

J \N,n) = (n-N/2)\N,n). (2.3) 

The operators J± will often be referred to as "ladder operators". Note that the action of J_ and J + 
on the end point vectors is J-\N, 0) =0 and J+\N, N) = 0. 
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It is convenient to introduce the number operator J\f, which is such that J\f\N,n) = n\N,n); it 
is easily seen that this operator is related to Jo by the formula N = Jq + N/2. The most general 
action of any powers of the ladder operators J± on the basis vectors is expressible in terms of the 
Pochhammer symbol. Indeed, one finds 



4\N, n) = J ^^^^ IN, n + k) = ^/(-l)*(n + l) fc (n - N) k \N, n + k), (2.4) 



n\{N -n + k)\ 



J " |iV ' n) = ^ (n-k)\(N-n)l lN,n ~ k) = ^-V k (- n ^ N ~ n + 1 )^ n ~ k )~ ( 2 - 5i 

These formulas are obtained by applying (|2.2p on the basis vector and by noting that \0d^)l 
(n + l) k and that = (-l) k (- n ) k [9]. 

2.1.2 Coherent states 

Let us introduce the su(2)-coherent states \N,rj) defined as follows 



I 1 N /AJ\V2 

where rj is a complex number. The action of the ladder operators on a specific coherent state \N, rj) 
can be computed directly to find: 

J+\N,7]) =7 1 - 1 N\N,ri), J-\N,rj) = 7](N -M)\N,7]). (2.7) 

These relations can be generalized to arbitrary powers of the ladder operators; one writes M in 
terms of Jq and uses commutation relations (12. lf> to obtain: 



J k + \N,7 1 ) = (-l) k i 1 - k (-M)k\N,7 1 ), 

J k \N, v ) = (-l) k rj k (Af-N) k \N,r,). (2.8) 
Note that the formulas involving the Pochhammer symbols are to be treated formally. 

2.1.3 Contraction of su(2) to rji 

The contraction of su(2) is the limiting procedure by which su(2) reduces to the Heisenberg-Weyl 
algebra. The Heisenberg algebra is generated by the creation-annihilation operators a', a and the 
identity operator. This algebra is denoted hi and defined by the commutation relations 

[a,a f ] = it), and [a, ill] = [a^iD] = 0, (2.9) 

where it) stands for the identity operator. The contraction corresponds to taking the limit as 
N — > oo; in order for this limit to be well defined, the ladder operators must be rescaled by a factor 
yN. The contracted ladder operators are denoted by 

j£°°> = lim i 



jV->00 y/N 
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In this limit, the action of the ladder operators on the basis vectors \N,n) tends to the action 
of the creation-annihilation operators in the irreducible representation of the rji algebra, which is 
infinite-dimensional. Indeed, the following formulas are easily derived: 

Vn + l\n + l), 

\fn \n — 1). 

This limit shall be used to establish the correspondence with studies associated to the Heisenberg- 
Weyl algebra; see |14j . 

2.2 Operators and their matrix elements 

The operators S, Q and their matrix elements, which will be the central objects of study, are now 
defined. 

2.2.1 S, Q, their matrix elements and biorthogonality 

Let a and b be complex parameters; we define 



with M 6 N. These operators can be represented by (N + 1) x (N + 1) matrices; their matrix 
elements are defined as 

V> M := (k,N\S\N,n), 
<p ki n := (k,N\Q\N,n). 

The two operators S, Q are obviously invertible, with inverses given by S- 1 = e- bJ -e- aJ + and 
Q^ 1 = e~ bJ - e~ aJ+ . We denote by Xn,k = (n,N\S~ 1 \N ) k) the matrix elements of <S -1 ; this leads 
to the following biorthogonality relation: 

N N 

J2xm,ki>k,n = ^(m,N\S- 1 \N,k)(k,N\S\N,n) = (m, N^ 1 S\N,n) = S nm , 

k=0 k=0 

where we have used the identity X^feLo 1-^' k)(k, N\ = it), which follows directly from the orthonor- 
mality of the basis {\N, k)}^ =Q . A similar relation can be written for the matrix elements of Q _1 . 
If one defines <; n ^ = (n, A^|Q _1 |A^, k), then 

N 
fc=0 



lim —p= 
lim , 



\N,n) 



\N,n) 
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2.2.2 Recurrence relations and polynomial solutions 

The algebraic nature of the operators S and Q leads to recurrence relations satisfied by the matrix 
elements ipk,n an d fk,n- Let us start by showing how the recurrence relation for the ipk,n s arises. 
The first step is to observe that 

(A; — N/2)ip kjn = (k,N\J S\N,n) = (k, NlSS^JoSlN, n). (2.10) 

The quantity S~ 1 JqS is computed using the Baker-Campbell-Hausdorff formula: 

e x Ye -x = Y + [XjY] + _L [X) [Xj Y]] + ±_ [Xj [Xj [Xt Y]]] + 

Applying this formula and the relations of appendix A, one readily finds 

S- 1 J S = J - 2bJ 2 _ + 2a[J+ + 26(1 + 2J ) J_ - 46 2 .fi] 2 . 
Substituting this result into (|2.10p yields for ipkn t ne recurrence relation 

2a v 7 (n + l) 2 {n - N) 2 ipk,n+2 = {k-n) tp k ,n + 2b\J (-n) 2 (N - n + 1) 2 ipk,n-2 

3 

—n)2t{N — n + l)2t 1pk,n-2ti 

t=0 

where the coefficients £t(n, N) are given by: 

£ (n,N) := 2[2n- N][2n 2 - 2nN - N+l), 
£l(n, N) := 4[6ra 2 - 6n(N + 2) + iV 2 + 5iV + 9], 

b(n,N) := 16(2n - N - 4), (2.11) 
Un,N) := 16. 

The recurrence relation for the matrix elements tpk,n is n °t OI the form (|1.1|) . However, the jumps 
on the index n are all even and it is clear from (I2.3P that any ipk,n with different index parity will 
be zero. Thus, setting n = 2j + q, k = 2t + q with q = 0, 1 and introducing the monic polynomial 

ipf,^ = ijjk,q y/x ^-qy.ni ) A^\£), the recurrence relation becomes 

A%{t) = (l- j)Af{l) + c[(-n) 2 (N -n + l) 2 ]^-i W 

3 

+ cJ2(-c)%(n, N) [(-n)at(JV - n + l) 2t ] ^- t (4 (2-12) 
t=o 

with c = a&. This relation is precisely of the form (|1.1P ; we thus conclude that the matrix elements 
ipk,n are given in terms of two families of cf-orthogonal polynomials with d = 3 corresponding to 
q = and g = 1 . These two families are fully characterized in the next section. 

The recurrence relation for the matrix elements (fikn can be obtained in the same fashion. It is 
first noted that 

(k-N/2)<p k , n = {k,N\J Q\N,n) = (k, NlQQ' 1 J Q\N, n). 
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The computation of Q 1 JqQ proceeds along the same lines. The result is 

Q- l J G Q = Jo + aJ + - MbJ M + abM(M - 1 + 2J )J M - 1 - ab 2 M 2 J 2 _ M - 1 . (2.13) 



Once again, introducing the monic polynomial ipk, n = Vifc.o ^ a n \j ^nw. J B n (k), this becomes 

B n+1 (k) = (k — n)B n (k)+f( M (n,N)B n „ M (k) 

+ f( M -i(n, N) B n+l - M + f <2M-i(n, N) B n +i-2M, (2.14) 

with / = a M b and where the coefficients £(n, N) are: 



Cm(u,N) 
( M -i(n,N) 
(2M-i(n,N) 



(-l) M M(-n) M (N-n+l) Mt 



(-l) M M(-n) M -i(N -n + l) M -l(2n -M-N + l), 

)2M-1- 



.lfM-l M 2 { _ n)2M _ i{N _ n + l) 



Therefore, the polynomials B n (k) are ci-orthogonal with d = 2M — 1. 
2.2.3 Contractions and the hi algebra 

It is relevant at this point to recall related results obtained in connection with the Heisenberg algebra 



f)X- These are to be compared with the contractions of the polynomials Aj (£) and B n (k) obtained 



in the next sections. 

In |15j . two of us investigated the matrix elements 



V ,M = (A; | e 6(a t ) 2 e- 2 |n>, 

where a and a) are the generators of hi and the vectors \n) are the basis vectors of its irreducible 
representation. It was shown that these matrix elements are given in terms of two series of Meixner 
polynomials with different parameters. Indeed, with n = 2j + q and k = 21 + q, one fmda3 

il>$ <x Mj(£; 1/2 + q,z). 

It is clear that this operator corresponds to the contraction of the operator S previously defined. 
Thus, the polynomials A^\k) are expected to tend to the Meixner polynomials in the limit as 
N — > oo and can be interpreted as a d-orthogonal finite "deformation" of Meixner polynomials. 
In |14j . we investigated the matrix elements 



e n) 



These matrix elements were found to be d-Charlier polynomials with d = M. These matrix elements 
correspond to the contraction of the matrix elements of the operator Q just defined. Consequently, 
it is expected that the contraction of the operator Q will lead to these d-Charlier. Moreover, the 
matrix elements <p k n are expected to yield the standard Krawtchouk polynomials K n ( X ; P , N) when! 
M = 1; the cases M/l correspond therefore to some d-Krawtchouk polynomials. 



See appendix B for definition and properties of Meixner polynomials M n (x; /?, z). 
2 See appendix B for definition and properties of Krawtchouk polynomials K n (x;p, N). 
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3 Characterization of the Aj(£) family 

We shall now completely characterize the family of d-orthogonal polynomials arising from the matrix 
elements of the operator S = e aJ + e bJ - ; these matrix elements have already been shown to satisfy the 
recurrence relation (|2.12p . The general properties are computed first and contractions are studied 
thereafter. 



3.1 Properties 

3.1.1 Explicit matrix elements 

We first look for the explicit expression of the matrix elements ipk,n = {k, N\S\N, n). This expression 
is obtained by setting n = 2j + q and k = 21 + q, expanding the exponentials in series, using the 
actions ([2~4")l and ([23]) and recalling the identity (a) 2n = 2 2n (f) { g ^ k ) n as well as (2n + q)\ = 
2 2n q\ nl (q + l/2) n . Extracting the factor 



_ (/ l(N-q)\k\ 

and pulling out the normalization factor ensuring that the polynomials are monic (a) - -' 
yields 



(JV-n)! 

(N-q)\n\ 



( (")~ J W A^ ™L )A ( f\£;c,N)^ k!q , (3.1) 



(N-q)\n\ 



where we have set 



3 1 ' ' ' ' j\ (N-n)\ 2 3 



-3 ~t J_ 
H +1/2 ^ 2z^±l'i 6c 



(3.2) 



2 2 

with c = ab. It is understood that if n and k have different parities, the matrix element ipk n is zero. 

We thus have an explicit representation of the polynomials Aj 1 \k;c,N) in terms of generalized 
hypergeometric functions. These polynomials satisfy the recurrence relation (|2,10p . Moreover, 
the fact that (k; c, N) is expressed as a indicates that these polynomials belong to the 
classification proposed in [3j. Indeed, it is clear that the singleton {—£} is 1-separable; consequently, 
the polynomials are of the form (ll.2p with s = 1. Thus, the polynomials A^\k; c, N) are examples of 
d-orthogonal polynomials corresponding to the case 1 of the classification with s = 1 and q = d = 3. 

3.1.2 Explicit inverse matrix elements 

The matrix elements Xn,k of the inverse operator S~ l = e~ bJ -e~ aJ + can also be evaluated explicitly; 
they can be computed either by directly expanding the exponentials or by inspection. Indeed, one 
finds the matrix elements of the inverse to be 

Xn,k = i>N-k,N-m (3-3) 

with * denoting the substitutions a —> —a and b —> —b. If N is an even number of the form 
N = 2p + 2q, the inverse is given by: 



Xn ' k { ' (p - ey.y k\(N - n)\ * ,c ' j ' 
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with n = 2j + q and k = 21 + q. If N is an odd integer, the cases q = and q = 1 have to be treated 
separately; since no further difficulty arise and Xn,k is expressed directly in terms of VVi,fc> we omit 
the details. 



3.1.3 Biorthogonality relation 

As pointed out in section 2.2.1, the matrix elements of the inverse 5 _1 provide the polynomials 
entering in the biorthogonality relations of the Aj (£; c, N) family. In view of fj3-3|) , this biorthog- 
onality relation reads 

N 

^2 ^k,n1pN-k,N-m = <W (3.5) 
k=0 

If N = 2p + 2q, the biorthogonality relation involves two Aj (£; c, N) with the same q; in the case 
where N = 2p + 1, the relation will involve two polynomials with different values of q. 
Firstly, set N = 2p + 2q and m = 2j' + q; the biorthogonality relation (I3.5P becomes 

1=0 

with the weight w e (p) = gr^Qi ■ 

Secondly, set N = 2p + 1 and m = 2j' + q; the biorthogonality relation is then of the form 

j2Mp)Af\£;c,N)A^ f (p-£;c,N) = (-l)^y. 
£=o 

A similar relation can be obtained with the q = and q = 1 polynomials in a different order. Thus, 
the two classes of polynomials corresponding to the values q = and q = 1 of the A^ 9 '' (£; c, iV) 
family are interlaced when the degree of the representation is odd and independent if the degree is 
even. 



3.1.4 Generating function 

We now derive the generating function for the A^ (£; c, N) polynomials. This derivation is based 
on the action of the operator S on coherent states \N,rj). Consider the following function: 

1 ^ /a-i(N-q)\ j2 fN\ 1/2 

n=0 

Upon substitution of the expression for the matrix elements, one finds 

N 



G(k-r ] ) = Y J Af{£;c,N) 



n! 

n=0 



setting G(k;rf) as a generating function for the polynomials A^ q \£;c,N). This generating function 
G(k;rj) is in fact the matrix element of S with respect to the coherent state \N, rj). Indeed, with 
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definition (12,61). it follows that 



N 



G(k; V ) = J_ >" g (^ g) ! Y^{k,N\S\N,n)(n,N\N,rj) 



Yli > q n=0 



1 a-i(N-q)\ ., - rl _ IAT . 

(fc, N\S\N,ri), 



i>k, q V iV! 

where the normalization factor was omitted. The matrix element (k, N\S\ N, rj) can be evaluated 
by using the spectral decomposition. One first decomposes the matrix element as: 

N 

(k,N\S\N,ri) = ^2(k,N\e aJ +\N,n)(i^,N\e bJ -\N,ri). 
The action of the ladder operators J± on the coherent states can then be used to obtain: 

L(tf-£)/2j 2m+q / ln/2 \ 

G(^;,) = ( e -/ 2 )^(_^)- £ ( ( ^| g) , (-^^- g ^j, (3.7) 
where 7J n is the standard Hermite polynomial 
3.1.5 Difference equation 

In a fashion dual to the way the recurrence relation was obtained, the difference equation satisfied 
by the matrix elements ipk,n can be derived with the help of the Baker-Campbell-Hausdorff formula 
and the formulas from appendix A. First, observe that 

(n- N/2)%b k , n = {k,N\SJ \N,n} = (k, N\SJoS~ 1 S\N, n). 
The operator SJ$S~ l must be evaluated here. Using the formulas given in the appendix, one obtains 

SJqS- 1 = Jo - 2aj\ + 26[J_ + 2aJ + (l + 2J ) - 4a 2 J%] 2 . 
Substituting the result, one finds 

(n - fc)^fc,n = 2by/ (k + l) 2 (k - N) 2 ^k+2,n ~ 2(L^ {-k) 2 (N - k + 1) 2 lp k -2,n 
3 

- 2ab^2(-a) t Uk,N) v / (-k) 2t (N - k + l) 2t ^_ 2t , n , 

where the coefficients £t(k,N) are those found in (|2,lip . The difference equation for the matrix 
elements ijj n k can straightforwardly be turned into a difference equation for the family of polynomials 

Af\k;c,N). Indeed, we have 

(j - £)2f {£; c,N)=c ^(q; N) If {£ + 1; c, N) - £A) q) (£ - 1; c, JV) 

3 

- c^2(-£) t Ct(2£ + g, AT) Af(£ - t; c, N), (3.8) 

t=0 



3 See appendix B for definition of the Hermite polynomials. 
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with 

Vt(N) = j^jW + q+ 1)2(21 + q- N) 2 . 

This difference equation can be written as an eigenvalue problem; denoting V/(x) = f(x) — f(x — l) 
Af(x) = f(x + 1) — f(x) and using the well-known identities 



f(x + t) = E(*) A W 

= E(-irQ)v 

on — n V / 



10=0 

the recurrence relation can be written as 



13 t /A 

H®=£V + c n e (q; iV) E A™ - ^("^(^ + ?, ^) E ^ ( w ) ' 



3.1.6 Shift operators 

The shift operators can also be derived from the Baker-Campbell-Hausdorff formula and the for- 
mulas from appendix A. We have 



V / (-n) 2 (A r -n + l) 2 Vfe,„-2 = {k,N\SJ 2 \N,n) = (k, N\SJ 2 S~ l S\N, n). 
One then easily obtains 

SJ 2 S- 1 = [J_ + 2aJ + (l + 2J ) - 4a 2 j|] 2 . 
Upon substitution in the previous equation, one gets 



{-n) 2 (N - n + l) 2 ^ k ,n-2 = V( k + 1M* ~ N) 2 ^k+2,n 

3 

- aJ2(- a )%( k > N W(-kht(N -k + l)2t V*-2t,n. (3-9) 



t=0 



where the coefficients N) are those given in (I2.1ip , Using the identities involving the difference 
operators introduced above, we obtain 

F^Af{i- c, N) = (-2j - q) 2 (N - 2j - q + 1) 2 A^i; c, N), 

with 

13 t /A 

= n e (q-, n) e a" - E(-a&(*, ^) E ( v " 

w=0 t=Q «j=0 ^ ' 

The computation of the backward shift operator for the Aj (£; c, N) is much more involved; the 
operator contains terms up to order 18. The explicit expression will not be given. 
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3.1.7 d-Orthogonality functionals 

From the recurrence relation (|2,10p computed in section 2, we concluded that the polynomials 
Aj(£;c,N) are (i-orthogonal with d = 3. We now provide the functionals with respect to which the 
Aj(£;c,N) polynomials are orthogonal. 

The initial step is to compute the recurrence relation satisfied by the matrix elements Xn,k of 
the inverse operator. First, we note that 

(k — N/2) X n,k = {^NIS^MN^) = {n^NlS^JoSS-^N,^. 

Proceeding along similar lines as before, we get 

S- 1 J S = J - 26J 2 + 2a[J+ + 26(1 + 2J ) J_ - 46 2 J 3 ] 2 . 

We then find the recurrence relation to be 

(k - n) Xn,k = 2aV '{-n) 2 (N - n + 1) 2 Xn-2,k - 2by/ (n + l) 2 (n - N) 2 Xn+2,k 

3 

-2a6^(-6)Vt(n,iV) v / (n + l) 2t(n - N) 2t Xn+2t,k, (3.10) 
t=o 

where the coefficients are: 

<7 (n, iV) := 2(2n - N)(l + 2n(n - N) - N), 
a x (n, 2V) := 4(6n 2 - 6n(iV + 2) + iV 2 - 7iV + 9), 
a 2 (n,N) := 16(2n - JV + 4), 
a 3 (n,iV) := 16. 

From this recurrence relation, the following proposition can be stated. 
Proposition 1. We /icwe 

2 
i=0 

where Sj(£) = X2i+(j,2£+q anc ^ 1^ is a polynomial in the variable £. The degree of the polynomial 
is determined as follows. Suppose j = 3j + 5 with 5 = 0, 1,2; then, if i < 5, the degree of the 
polynomial is equal to 7, otherwise it is equal to 7 — 1. 

Proof This statement is obtained directly from the recurrence relation. Setting n = 2j + q = 
yields the element X2(3)+q,k in terms of the elements X2(2)+ q ,k, X2(i)+ q ,k, Xg,k with degree 1 in £ in 
the case of Xo+q,k and for the others. The proposition then follows by induction. □ 

With this statement, we define the following linear functionals 

L(JV-q)/2j / — 

for i = 0, 1,2. From the biorthogonality relation, it follows that 

C^Af{i;c,N)][=^ ii 3>*T + i + 1 -> (3 . n) 
1 L J v 7J 1/0, if j = 37 + 2. V ; 
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With these relations, the d-orthogonality of the family of polynomials Aj(£; c, N) is manifest. 

The algebraic setting at hand has thus allowed to completely characterize the family of polyno- 
mials Aj (£; c, N). The Baker-Campbell-Hausdorff relation and the formulas from appendix A were 
used to derive the recurrence relation, the difference equations as well as the forward shift operator. 
Moreover, the action of the operators J± on the basis vectors \N,n) of the (N + l)-dimensional 
irreducible representation of su(2) was used to compute explicitly the expression for the polynomials 
Aj(£\c,N). As it turns out, these polynomials are d-orthogonal with d = 3. They have a simple 
representation as 2-^3 generalized hypergeometric function and fall into the classification of Ben 
Cheikh et al. due to the 1-separability of the singleton set {£}. The inverse elements were found us- 
ing the symmetry of the vectors \N, n) with respect to the action of the operators J±; those inverse 
elements proved useful in obtaining the three orthogonality functionals for the Aj(£; c, N)'s. In the 
next subsection, we will look at the contractions of the Aj {£; c, N) and their structural formulas. 

3.2 Contractions 

We now turn to the evaluation of the contraction of the Aj(£; c, N). We recall that the procedure 
of contraction, explained in section 2, corresponds to taking the limit as N — > 00 after renormalizing 
the operators J± by a factor of yN. 

3.2.1 Contraction of the recurrence relation 

Let us first apply the contraction to the recurrence relation. The renormalization of the su(2) 
generators is implemented by substituting c by c/N 2 . We denote 

Jim Af(£;c,N) = M J (£), 

provided that this limit exists. After straightforward manipulations, one finds that after the renor- 
malization, the terms in Ay}_ 3 (£; c, N) and A^} 2 (£; c, N) are of order 0(N~ 2 ) and 0(N^ 1 ), respec- 
tively. Consequently, they tend to zero in the limit as N — > 00 and the recurrence relation becomes 

£ Mj(£) = M j+ i(£) + {j - c [2 + 4(2 j + q)]]^) 
- [c(2j + q)(2j + q - 1)[1 - 4c]} Mj-i{£). 

The initial 5-term recurrence relation contracts to a 3-term recurrence relation. To identify to which 
orthogonal polynomial this relation corresponds, we set -4c = j^. The recurrence becomes 

£ Mj {£) = M j+1 (£) + -^{j + d(j + q + 1/2) }m^) 

+ (T^{ (j + q/2){j + q/2 " V^j^-xW. (3.12) 

We recognize in (|3.12p the normalized recurrence relation of the Meixner polynomials Mj(£; (3,d) 
with P = q + 1/2, for q = 0,1. Thus, as expected, we have 

lira A) q) (£; c, N) = Mj q + 1/2, -?-) . 

N^oo J \ C — 4 J 
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3.2.2 Contraction of the matrix elements 



It is relevant to examine directly the contraction of the explicit formulas obtained for the matrix 
elements ip n k- To take this limit, one must first expand ip n k by writing the generalized hyper- 
geometric function as a truncated sum. Then, using the same renormalization as in the previous 
contraction, one straightforwardly obtains 



lim ip n k 

N^oo 



-j-l 1_ 
( 9 + l/2)'4c 



(3.13) 



Here, the explicit expression of the Meixner polynomials in terms of Gauss hypergeometric functions 
is recovered. The contraction at this level exhibits clearly the relationship between the parameters 
of the polynomial family A, (£;c, N) and the Meixner polynomials. We have (3 = q + 1/2 and 
1 — \ = i-. This is the result found in [15j@ 



3.2.3 Contraction of the generating function and coherent states 

The contraction limit of the generating function G(k; rf) can also be taken. However, it is easy to see 
that in the limit N — > oo, the coherent states \N, rj) as given by (|2,6p are ill-defined. Indeed, under 
the contraction of the su(2) algebra, the radius of the Bloch sphere, on which the coherent states 
are defined, must also be taken to infinity. Therefore, a well-defined contraction of the coherent 
states requires to take the limit N — > oo with the renormalization r] — > rj/yN |11| . Note that with 
this renormalization, the coherent states become eigenstates of the operator J-/yN in the limit 
N -> oo. 

Performing the same renormalization, the contraction of the generating function yields that of 
the normalized Meixner polynomials 

KmG(*;»7) =jr,M j (£,q+ 1/2, c)^^ = rf 1 F 1 

°° n=0 

3.2.4 Contraction of inverse matrix elements and biorthogonality 

In section 3.1, the symmetry of the irreducible representation of su(2) was used to obtain that 
Xn,k = ^N-k N-n- This result is due to the fact that in this representation, there exists a vector 
\N, 0) which is annihilated by J_ and a vector |JV, N) which is annihilated by the operator J + . In the 
contraction limit, the symmetry of the representation is not preserved. Indeed, the vectors of the f)i 
irreducible representation are also labeled by a positive integer n, but this integer is unbounded from 
above; more precisely, there exits no basis vector such that a)\n) = 0, but the relation a\0) = still 
holds. Consequently, all relations involving the inverse matrix elements Xn,k must be recalculated 
with the contraction applied directly to the operators. We stress that this loss of symmetry under 
the limit N — > oo is the reason for the drastic change of behavior of the polynomials Aj (£; c, N) 
and for the fact that, in particular, the d-orthogonality reduces to the standard orthogonality. 

To find the expression for the inverse matrix elements in the contraction limit, one can simply 
calculate, directly from the operators, the matrix elements of S* -1 = e~ ba2 e~ a ( at ) 2 written as Xn,k = 

4 To recover the full result, one has to take a — > —ZJ, b — > lj and eliminate the diagonal term in expansion of |15| . 



-I rf 

g+1/2' 4^ 
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(n\S~ 1 \k) or take the contraction of the recurrence relation. Meixner polynomials with changes in 
the arguments are found. 

Other limits involving objects such as the difference equation, the forward shift operators can 
be taken in the same way, yielding their counterparts for the Meixner polynomial. 

4 Characterization of the B n (k) family 

In this section, we fully characterize the d-orthogonal polynomials, with d = 2M — 1, in terms of 
which the matrix elements of the operator Q = e aJ+ e bJ - are expressed; these polynomials have 
already been shown to obey the recurrence relation (I2.14p . 



4.1 Properties 

4.1.1 Explicit expression for the matrix elements of Q 

The matrix elements (pk, n = {k, N\Q\N, n) can be computed explicitly by expanding the exponentials 
in series and using the actions (|2.4p and (|2.5p . To express the matrix elements in terms of generalized 
hypergeometric functions, one needs the identities 



(a) 



Mn 



M 



Mn TT / ° + P 



n 

0=0 



M 



M-l 



and (Mn + q)\ = M Mn q\ JJ 



/3=0 



q + (3 + l 
M 



Setting n = Mj + q with q £ {0, . . . , M — 1}, one finds 
a k -ib> 



ak - q V ^ (N-q)l(N-q)\ 
j\{k-q)\q\ ' "I/ (^-fc)!(JV-n)! 



\ M I -2M 



~j {«m} -1 
{/3 m } {7m}' (Ma) M b 



with {a m } = (q — k + m)/M, {/3 m } = (q + m + 1)/M with q + m + 1 = M excluded from the 
sequence and {"f m } = (q — N + m)/M, with m running from to M — 1. 

To obtain the exact expression for the B n (k) polynomials, one must pull out the "ground state" 



N 



1/2 



and the normalization factor a n y ^,^ ! from the expression of the matrix elements. The final 
expression reads, with / = a M b: 

Vk,n = B n (k;f,N) 

with 



Wnl 



B n (k;j,J\) = — t— t 1+M-T2M-1 

•1/7 



{a m }. 



-1 



{Pm} {7m}' M M f 



(4.1) 



j\q\ (N — n)\ 

The recurrence relation (I2.14|) indicates that these polynomials are d-orthogonal with d = 2M — 1, 
but it is clear that the set {a m } is s-separable only for M = 1, it is 1-separable in this case. 
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However, with this value of M, the polynomials are expressed as Gauss hypergeometric functions 
and are simply related to the Krawtchouk polynomials. For M = 1, the matrix elements are 

1/2 / . A 1/2 



<p k ,n = a k b n (^ n J (J K n (k;p,N), 



with p = —ab. For any other value M, the polynomials are (f-orthogonal extensions of the 
Krawtchouk ones, but they fall outside the classification of [3]. 

4.1.2 Matrix elements of Q 1 

The matrix elements of the inverse operator ^ n k = (n, A^Q" 1 ] N, k) with Q~ l = e~ bJ - e~ aJ+ can 
also be computed directly. Just as for the matrix elements ifj n ^ of the S operator, the matrix 
elements <; n ^ of the inverse operator Q~ l possess a reflection symmetry; we indeed find 

?n,fc = tPN-k,N-n> ( 4 - 2 ) 

where * denotes the replacements a — > —a and b — > —b. In terms of the matrix elements (fk, n , this 
biorthogonality relation reads 

N 
fc=0 

Because of the asymmetric form of the operator, the explicit expression for the matrix elements 
heavily depends on the behavior of k and N modulo M. An exact expression would thus comprise 
M 2 cases of residues. 

For definiteness, we shall set M = 2 whenever general expressions cannot be found in closed- 
form. Note that this case is relevant in the contraction limit because such powers of the creation- 
annihilation operators appear in the oscillator algebra scf^ |15| . 

4.1.3 Biorthogonality relations 

The biorthogonality of the B n (k; /, N) polynomials can be written for any value of M. Indeed, one 
has 

N 

^ w k B n (k; f, N)B N ^ m (N - k; /', N) = (-l) n 5 nm , (4.3) 

with Wk = fc| (jv-fc)! an< ^ f = (~^) M+1 f ■ Note that when M = 1, this is not exactly the orthog- 
onality relation of the Krawtchouk polynomials. To obtain the orthogonality of the Krawtchouk 
polynomials from this equation, one must remove the normalization factor and use Pfaff 's transfor- 
mation. This transformation is only available for 2-F1 hypergeometric functions, thus, M = 1 is the 
only case for which the biorthogonality relation degenerates into the standard orthogonality. 

4.1.4 Generating function 

The generating function for the B n (k; /, N) polynomials is obtained as in the previous section; it 
can be derived explicitly for any value of M. Define 



N 



aV /2 



G(k; V ) = <p k , n rf. (4.4) 
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Substituting the expression for the matrix elements, one gets 

JV 



G{k;7 ] ) = Y,B n (k-f,N) {ri/i ' 



n=0 



n 



yielding a generating function for the polynomials B n (k; f, N). Once again, this generating function 
is expressed as the overlap between the vector \N, k) and Q\N,rj). Indeed, it follows from the 
definition of the coherent states (12. 6ft that 

G(k; V ) = —(k,N\Q\N, V ). 
Using the action of the ladder operators (12. 7|) and (12.8|) . one easily finds 

G(k; V ) = mF [{S m };(-l) M (Mr,) M b] , (4.5) 

with {S m } = ^~jvf +m with m = 0, . . . , M — 1. After a slight adjustment in the parameters and use 
of the identity \Fq(— a; b) = (1 — b) a , the generating function for the Krawtchouk (see Appendix B) 
polynomials can be recovered for M = 1. 

4.1.5 Difference equation 

The algebraic setting can be used to derive the difference equation satisfied by the matrix elements 
(Pn t k- Again, one writes 

(n-N/2) VKn = (k,N\QJ \N,n) = {k, NlQJoQ^Q^n). 

From the Baker-Campbell-Hausdorff relation, it follows that 

QJoQ^ 1 = J - aJ + + Mb[J- + 2aJ - a 2 J+] M . 

Setting M = 2, the difference equation for the matrix elements is found to be 

(n - k) ip k , n = 2by/(k + l) 2 (k - N) 2 <p k +2,n + 4ab y/(k + 1)(N - k) Ci <Pk+i, n 



+ 2 a 2 b Co <Pk,n ~ 4a 3 by/k(N - k + 1) Ci Vh-i,n 

- a^/k(N - k + 1) v9 fc _ ljn + 2a 4 b v / (-k) 2 (N - k + 1) 2 ^ fe _2,„, (4.6) 

with the coefficients Ci : 

C = 6fe 2 -6/c7V + A^(A + l), 
Ci = 2k - N - 1. 

For the B n (k) polynomials, this equation becomes 

nB n (k; /, AO = (2f)(k - N) 2 B n (k + 2; /, N) + (4 / Ci)(AT - fc)/^ + 1; /, N) 
+ (k + 2f ( )B n (k; f, N) - k(4f Ci + l)B n (k - 1; /, N) 
+ (2f)(-k) 2 B n (k-2;f,N). (4.7) 
Using the same identities as before, this difference equation can be written as an eigenvalue equation. 
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4.1.6 Orthogonality functionals 

The d-orthogonality functionals can be computed for the polynomials B n (k;c,N). First we have: 

(k-N/2)^ k = {n,N\QJ \N,k) = (n, N\QJ Q Q~ l Q\N, k). 

For M = 2, the formula (pTl3]) yields 

Q- 1 J Q = J - 26J 2 + aJ + + 2o6(l + 2J ) J_ - 4a6 2 j£. 

Substitution in the first equation gives a recurrence relation for the matrix elements of the inverse 
operator Q _1 ; this relation is 



(k - n) ^ k = a^n(N -n + l) <^_ lifc + 2a6(2n + 1 - N)^{n + 1)(N - n) ?n+ljfc 

- 26^(n + l) 2 (n-iV) 2 ?n+2)fc - 4a6 2 y/-(n + l) 3 (n - iV> fe+3 , t . (4.8) 

The form of this recurrence relation suggests the following proposition: 
Proposition 2. We have 

2 

, n>k = J2 Y t ] mi(k), (4.9) 
i=0 

with Hj(fc) = q ^ and a polynomial in k. If n = 37 + 5 with S = 0, 1, 2; i/ten i/ie degree of the 
polynomial is 7 u>/ien i ^ 5 and 7 — 1 otherwise. 

Proof The proof follows from the recurrence relation and by induction. □ 
This suggests the definition of the linear functional 

- //V\ 1/2 

Mrf(x) = Y / a X ( ^(x)/(x). (4.10) 

It thus follows that 

M,^Mt,c,N))H + i + 1 (4.11) 
I 7^ if n = 37 + i 

General orthogonality functionals could be defined in full generality for the operator Indeed, 
one could emulate the procedure to build 2M — 1 orthogonality functionals. 

4.2 Contractions 

The contractions of the polynomials B n (k) are expected to yield the d-Charlier polynomials con- 
sidered in |14j . The loss of symmetry also occurs in the contractions of the polynomials B n (k); 
consequently, contractions of the orthogonality functionals and ladder operators cannot be taken 
directly. As those cases were treated in full generality in |14j . we simply show how the matrix 
elements (pk n contract to the particular case of (i-Charlier polynomials that was considered. 
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4.2.1 Contraction of <pk 



The required renormalization is a — >■ a/y/N and b — > b/VN M . Once again, the limit is taken by 
first expanding the generalized hypergeometric function as a truncated sum and performing the 
indicated substitution. Upon simple transformations, the result is found to be 



lim tp k r 



l k -4V 



:Vk\n\ i +M F M - 



j\(k - q)lq\ 

which gives precisely the matrix elements obtained in [14 



-3 (~1) M+1 
{An} ~ ' a M b 



(4.12) 



5 Conclusion 

We studied the d-orthogonal polynomials related to the classical Lie algebra 5u(2). We showed 
the matrix elements of the operators S = e aJ +e bJ - and Q = e aJ+ e bJ - were given in terms of two 
families of polynomials Aj {£;c, N) and B n {k; c, N). Using the algebraic setting, we characterized 
these polynomials; their explicit expressions in terms of hypergeometric functions allowed to identify 
those that belong to the classification given in [3]. 

We also studied the contraction limit in which su(2) is sent to the Heisenberg algebra fj^. We 

showed that the Aj (£; c, N) tend to standard Meixner polynomials when N — > oo. The Aj (£; c, N) 
can therefore be seen as discrete d-orthogonal versions of the Meixner polynomials. In addition, it 
was shown that the polynomials B n (k;c,N) are some d-orthogonal Krawtchouk polynomials and 
that they converge to the (i-Charlier polynomials in the contraction limit. 

In |15j . the exponentials of linear and quadratic polynomials in the generators of the f)i algebra 
were unified to yield matrix orthogonal polynomials; these considerations were motivated by the 
link with the quantum harmonic oscillator. Similarly, considerations regarding the discrete finite 
quantum oscillator would suggest to convolute the corresponding matrix elements. While providing 
physically relevant amplitudes for the quantum finite oscillator, this study could lead to <i-orthogonal 
matrix polynomials, such as considered in |12j . We plan to report elsewhere on this. 



Acknowledgments 

A.Z. wishes to thank the CRM for its hospitality. V.X.G. has benefited from a scholarship from 
Fonds quebecois de la recherche sur la nature et les technologies (FQRNT). The research of L.V. 
is supported in part through a grant from the Natural Sciences and Engineering Research Council 
(NSERC) of Canada. 



Appendix A— Useful formulas for su(2) 

The relation 

J Q Jl = J^(Jo±n), 

holds and can be proven straightforwardly by induction on n. Using this identity and the relations 
(J±y = J=p as well as Jq = Jo , it follows that for Q(J±) denoting a polynomial in J±, one has 

[Q(J±),Jo] = tJ±Q'(J±), 
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where Q'{x) denotes the derivative with respect to x. The preceding formula and the Baker- 
Campbell-Hausdorff relation lead to the identity 

e w±)j oe -w±) = j tj±q'(j±). 

In addition, we have the relations 

[J+, JH] = 2nJ G ,p- 1 + n(n - 

[J_, J£] = -2nJ^ _1 J - n(n - 1) J£~\ 

which can also be proved by induction on n. With the help of the previous identities, one obtains 

[J+,Q(J-)\ = 2J Q'(J_) + J-Q"(J-)> 
[J_,Q(J+)] = -2Q'(J+)J - J + Q"(J+), 

From these formulas it follows that 

e Q(J -^ J+e- Q{J -^ = J+- 2J Q'(J_) - J_[Q"(J_) + Q'(J-) 2 ], 
e Q(J+) j_ e -Q(J+) = j + 2Q'(J + )J + J+[Q"(J+) - Q'(J+) 2 ]. 

Appendix B Meixner, Hermite and Krawtchouk polynomials 
Meixner polynomials 

The Meixner polynomials have the hypergeometric representation 



M n (x;(3,d) = 2 F 1 



They satisfy the normalized recurrence relation 



-n, —x 1 



xB n (z) = B n+ i(x) + ) / ; B n (g) + v . *V _ ; ^n-l(x), 



1-d 



(i - ay 



with 



M n {x-p,d) = -l-( d —±\ B n (x). 



W)n V d 



Hermite polynomials 

The Hermite polynomials have the hypergeometric representation 

n 1 — n 

H n (x) = (2x) n 2 F ' " 



2' 2 : -i 
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Krawtchouk polynomials 

The Krawtchouk polynomials have the hypergeometric representation 



K n (x;p,N) = 2 F 1 



—n, —x 1 



-N 



P 



Their orthogonality relation is 




£ )p x (l - P ) N - x K m (x;p, N)K n {x;p, N) 



(-l) n n! fl-p 
H^T \~P~ 




They have the generating function 




n 



For further details, see [8]. 
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